DISTRIBUTION  STATEMENT  A.  Approved  for  public  release;  distribution  is  unlimited. 


Surface  Wave  Dynamics  in  the  Coastal  Zone 

Gerbrant  Ph.  van  Vledder 

Department  of  Civil  engineering  and  Geosciences,  Delft  University  of  Technology 
P.O.  Box  2048,  2600  GA  Delft,  The  Netherlands 
t:+31  15  2783255,  f: +31-15  2784842,  e:  g.p.vanvledder@tudelft.nl 

& 

BMT  ARGOSS  bv,  PO.  Box  61,  8325  ZH,  Vollenhove,  The  Netherlands 
t  :+3 1-527-242299,  f:+3 1-527-2420 16;  e:  gerbrant.vanvledder@bmtargoss.com 

Leo  H.  Holthuijsen 

Department  of  Civil  engineering  and  Geosciences,  Delft  University  of  Technology 
P.O.  Box  2048,  2600  GA  Delft,  The  Netherlands 
t:  +31  15  2784803,  f:  +31-15  2784842,  e:  l.h.holthuijsen@tudelft.nl 

Kevin.  C.  Ewans 

Shell  International  Exploration  and  Production  B.V. 

Regional  Discipline  Lead  -  Metocean,  Sarawak  Shell  Bhd, 

50450  KUL-D WO,  Malaysia 
e:  Kevin.ewans@shell.com 

Award  Number;  N00014-12-1-0534 
http  ;//www. swan.tudelft.nl 


LONG-TERM  GOALS 

The  proposed  work  will  contribute  to  the  improvement  of  existing  third-generation  (3G)  wave  models 
as  well  as  to  the  development  of  the  next  generation  of  numerical  wave  modeling  capability.  The 
results  will  be  applicable  in  the  coastal  zone  from  deep  water  up  to  and  including  the  surf  zone.  Our 
efforts  will  focus  on  analyzing  high  quality  datasets  to  support  further  development  of  the  source  terms 
for  triad  interactions  (Snl3),  depth  induced  wave  breaking  (Sbrk)  and  bottom  friction  (Shot)  in  the 
near-shore  zone.  Another  point  of  interest  is  the  generation  of  bounded  long  waves  in  the  surf  zone. 

OBJECTIVES 

The  scientific  or  technological  objectives  of  this  project  are  to  understand  the  physical  processes  of  the 
evolution  of  wind  waves  in  the  coastal  zone  and  develop  accurate  parameterisations  of  these  processes 
for  application  in  numerical  wave  prediction  models. 

APPROACH 

The  proposed  work  is  subdivided  in  five  main  work  packages  (WP). 

1 )  Assembly  of  high  quality  data  set; 

2)  Analysis  of  spectral  evolution; 

3)  Development  of  a  source  term  for  wave  breaking  in  shallow  water; 
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4)  Development  of  a  source  term  for  triad  interactions. 

5)  Improvement  of  source  terms  for  bottom  friction. 

WORK  COMPLETED 

Collection  of  data  sets 

TU  Delft  and  BMT  Argoss  extended  their  collection  of  shallow  water  wave  data  suitable  for  the 
calibration  and  verification  of  newly  developed  source  terms  for  wave  breaking,  bottom  friction  and 
triad  interactions.  This  set  comprises  wave  flume  and  field  data.  An  overview  of  the  most  relevant  data 
set  collected  is  given  in  the  previous  annual  report.  Full  details  are  provided  in  the  manuscripts 
prepared  by  Salmon  and  Holthuijsen  (  2013,  and  Salmon  et  al.  (2013). 

Development  of  source  term  for  wave  breaking 

An  study  of  the  literature  on  depth-induced  wave  breaking  was  carried  out  resulting  in  a  comparison  of 
12  different  parameterizations  using  4  different  dissipation  models  (and  versions  thereof)  summarized 
in  Figure  1 . 


Figure  1.  Scatter  index  of  the  12  models  based  on  13  data  sets  containing  225  cases  of 
laboratory  observations  (*)  and  field  observations  (**).  The  highlight  colors  indicate  the  two 
ranges  of  performance  (green  and  blank)  and  four  ranges  of  scatter  index  (from  blank  to  orange). 

The  main  conclusion  of  this  comparison,  in  agreement  with  previous  studies  e.g.  Apotsos  et  al.  [2008], 
is  that  no  one  model  with  default  settings  provided  the  best  predictions  for  significant  wave  height. 
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Furthermore,  even  the  simplest  models  e.g.  Battjes  and  Janssen  [1978]  and  Rattanapitikon  [2007] 
performed,  on  average,  similarly  to  more  recent  and  complex  parameterizations  e.g.  Filipot  and 
Ardhuin  [2012]. 

This  disappointing  conclusion  led  to  the  P  —  kd  scaling  for  the  Battjes  and  Janssen  dissipation  model 
(Fig.  2)  which  provides  a  joint  dependency  on  both  local  slope  and  local  water  depth  normalised  by 
wave  length.  This  new  scaling  was  shown  to  resolve  the  problem  of  overestimating  wave  heights  in 
horizontal  laboratory  cases  e.g.  Jensen  [2002]  and  underestimating  wave  heights  in  (1-D  idealized) 
cases  of  locally  generated  waves  e.g.  Young  and  Babanin  [2006]  and  Bottema  and  van  Vledder  [2009]. 
However,  the  performance  of  the  p-kd  over  the  field  cases  showed  no  significant  improvement,  and 
in  fact  a  reduction,  in  the  skill  of  predicting  significant  wave  heights  (Fig.  3). 


Figure  2.  Calibrated  as  a  function  of  bottom  slope  n  =  tan  '  p  and  normalized  wave 
number  kd  and  ygj  =  0.73  for  reference. 
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Figure  3.  Scatter  index  for  the  verification  of  the  p-kd  scaling,  with  or  without  directional 

partitioning  for  cr*  =10°, 15°, 20°, 25°  and  30°  (results  with  a  fixed  scaling  =0.73  as  reference 

only).  Performance  over  laboratory  observations  (*)  and  field  observations  (**)  are  highlighted 
in  three  ranges  of  scatter  index  (from  blank  to  dark  gray).  The  Lakes  data  set  (used  in  the 
calibration)  is  not  included  in  computing  the  average  values  but  are  shown  to  demonstrate  the 
effect  due  to  directional  spreading  for  this  data  set.  Laboratory  cases  with  their  long-crested 
waves  are  unaffected  by  directional  partitioning  and  are  shown  with  a  cross-hatched  background 
where  directional  partitioning  is  used. 


It  was  speculated  that  for  field  cases,  the  analogy  of  Battjes  and  Janssen  [1978]  of  a  1-D  bore  should 
not  necessarily  hold  true,  particularly  for  short-crested  seas  topographically  induced  e.g.  by  a 
submerged  shoal  and  that  by  partitioning  the  wave  field  using  a  characteristic  directional  width,  the 
dissipation  due  to  depth-induced  breaking  could  be  better  represented.  This  effect  was  added 
heuristically  in  the  forms  of  Eq.  1  and  2: 
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With  the  number  of  wave  partitions  estimated  with  Kg=(jgl  (jg>\  where  a*g  represents  a 
eharaeteristic  direetional  width  over  whieh  the  1  -D  bore  analogy  can  be  assumed  to  be  appropriate. 
The  use  of  this  directional  partitioning  ( cr^  =  1 5°  )  with  the  P-kd  scaling  has  been  shown  to 
significantly  improve  the  prediction  for  significant  wave  height  in  horizontal  field  cases  (see  Fig.  3). 
Currently  the  effect  of  adding  wave  directionality  to  the  2  best  performing  models  for  field  cases 
(Battjes  &  Janssen  [1978]  with  y^j  -  0.73  and  Baldock  et  al.  [1998]  corrected  by  Janssen  and  Battjes 
[2007]  and  Alsina  and  Baldock  [2007]  with  the  scaling  of  Ruessink  et  al.  [2003];  see  Fig.  1)  is  being 
investigated.  Preliminary  results  show  that  these  2  parameterizations  and  the  P-kd  scaling  can  yield 

on  average  similar  or  better  results  in  the  field  cases  considered  if  15°  <  cr^  <  20° .  However,  the 
P  —  kd  scaling  still  on  average  performs  best. 


Development  of  source  term  for  shallow  water  triad  wave-wave  interactions 

An  inconsistency  was  identified  for  the  implementation  of  the  LTA  and  DCTA  models  in  SWAN 
arising  from  the  extension  to  2D  spectra  of  ID  source  terms  i.e.  (/)  to  5^,3  (/,^)  (required  for 
wave  models  dealing  with  2-D  energy/action  densities)  due  to  the  incorrect  assumption: 

\E{p,e)dex\E{p,e)de=\E{p,e)E{p,e)de 


This  appears  to  have  led  to  the  default  calibration  parameter  for  the  LTA  model  to  have  been  set  at  a 
low  value  of  0.05  (historically  0.5  and  0.1)  when  Eldeberky  [1996]  suggests  a  value  of  0(1). 

To  avoid  an  integral  over  all  directions  per  frequency  and  directional  bins,  which  would  substantially 
increase  computational  costs,  for  1-D  laboratory  cases  where  all  the  energy  are  located  in  only  a  few 
directional  bins,  the  source  term  for  the  2-D  spectrum  for  the  LTA  source  term  (sum  contribution)  can 
be  approximated  as: 


KlAf^^)-^LTAX 


f  ' 
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where  represents  a  calibration  parameter  of  0(1)  and  %  represents  the  remaining  terms  of  the 


original  expression  of  Eldeberky  [1996]  (his  Eq.  7.25).  A  similar  approach  was  used  in  the 
implementation  for  the  source  terms  for  the  SPB  (Becq-Girard  et  ah,  1999),  Toledo  and  Agnon  [2012] 
and  DCTA  model  (in  energy  conservative  form)  (Booij  et  ah,  2009).  These  implemented  source  terms 
were  then  calibrated  over  a  number  of  laboratory  cases  resulting  in  new  calibration  parameters  of  0.65, 
0.75,  0.35  and  0.50  for  the  LTA,  SPB,  energy  conservative  DCTA  and  Toledo  and  Agnon  models 
respectively.  In  addition,  the  original  coelficients  given  by  Becq  et  al.  were  re-assessed  by  setting 
- 1 .  The  models  were  then  verified  against  more  laboratory  observations  (Pig.  4). 


As  shown  in  Pig.  4,  by  using  bulk  parameters  alone,  it  is  not  possible  to  select  a  better  model  between 
the  calibrated  ETA,  DCTA  and  SPB  models.  However  both  the  DCTA  and  SPB  models  appear  to  give 
better  spectral  shapes  than  the  ETA.  Current  work  is  now  focused  on  using  a  performance  indicator 
based  on  the  transfer  of  energy  from  the  primary  peak  to  the  higher  harmonics  based  on  the  approach 
of  Becq  et  al.  [1999]  to  select  the  best  model  to  further  develop. 
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Model 

LTAo 

LTA 

SPB^ 

DCTAe 

SPB^ 

Beji-Battjes 

0.18 

0.18 

0.20 

0.18 

0.21 

Boers 

0.15 

0.14 

0.29 

0.11 

0.14 

LOWISH 

0.09 

0.10 

0.07 

0.09 

0.07 

Smith 

0.10 

0.10 

0.11 

0.11 

0.12 

Wallingford 

0.09 

0.08 

0.10 

0.11 

0.08 

Average 

0.12 

0.12 

0.15 

0.12 

0.12 

Figure  4.  Average  scatter  indices  of  and  ^  from  the  verification  over  laboratory 

observations  of  the  new  calibration  parameters.  The  current  default  for  SWAN,  LTAo,  is  added 
for  reference  The  model  of  Toledo  and  Agnon  are  not  included  due  to  its  undeveloped  state  and 
poor  performance  in  shallow  water  {kd<  1)  (e.g.  Fig.  5). 

^Represents  the  SPB  model  with  =  0.75  and  K  =  Q.95kp  -0.75  i.e.  using  the  coefficients 
suggested  by  Becq  et  al.  and  calibrating. 

^Represents  the  SPB  model  with  a^pg  =1.00  and  K  =  0.70/:^  -0.10  i.e.  re-calibrating  he 
coefficients  of  Becq  et  al.  (1999). 


Future  work  will  involve  verifying  an  improved  LTA,  DCTA  or  SPB  model  for  applicability  with 
third-generation  wave  models  in  field  conditions  and  reducing  the  computational  effort  required. 


Figure  5.  The  evolution  of  the  spectrum  (stations  1,  2,  5  and  7  in  the  upper  panel),  computed 
with  the  LTA  (Eldeberky,  1996),  DCTAe  (Booij  et  al.,  2009)  and  SPB^  (Becq-Girard  et  al.,  1999) 
and  Toledo  and  Agnon  [2012  Jtriad  source  terms  over  a  1:30  flatbed,  compared  with  the  spectra 
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observed  hy  Smith  [2004].  Note  that  the  vertical  scale  is  such  that  a  yt^'^-tail  would  appear  as  a 
horizontal  line. 

Long  wave  generation  in  the  coastal  zone 

The  SWASH  model  has  been  applied  to  simulate  long-wave  generation  at  the  FRF,  Duek,  NC.  The 
eomputational  results  were  eompared  with  observations  and  show  reasonably  agreement  with 
measured  data.  The  experiement  were  exeeuted  with  one  layer  in  the  SWASH  model.  We  are  now 
investigating  the  influenee  of  additional  layers  in  the  vertieal  on  the  predietion  of  long  waves  in  the 
eoastal  zone.  The  eomputations  were  exeeuted  at  Imperial  College  London,  in  elose  eollaboration  with 
Shell  and  Delft  University  of  Teehnology. 

RESULTS 

The  eolleetion  of  shallow  water  wave  flume  data  has  been  extended  and  rearranged  with  additional 
data  for  the  verifieation  of  new  parameterisations  of  shallow  water  proeesses. 

A  new  parameterisation  for  wave  breaking  in  shallow  water  over  sloping  and  flat  bottoms  was 
reformulated,  ineluding  a  dependenee  on  direetional  spreading  and  implemented  in  the  SWAN  model. 

A  new  parameterisation  for  triad  interaetions  in  shallow  is  being  developed.  As  a  first  step,  some 
ineonsisteneies  in  the  LTA  formulation  were  removed.  In  a  seeond  step  the  method  of  Beek  (1999)  has 
been  implemented  and  promising  results  were  obtained. 

The  SWASH  model  has  sueeesfully  been  applied  to  prediet  the  generation  of  long  waves  in  the  eoastal 
zone. 
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IMPACT/APPLICATIONS 


Economic  Development 

The  improvements  to  coastal  wave  prediction  models  as  developed  in  this  project,  will  contribute  to 
various  economic  activities  on  the  continental  shelf  and  the  coastal  zone,  such  as  fisheries,  shipping 
harbor  development  and  offshore  industry.  Further,  the  availability  of  improved  nearshore  wave 
prediction  models  will  benefit  coastal  and  ocean  engineering  companies  e.g.  in  the  design  and 
operation  of  offshore  and  coastal  structures,  and  the  development  of  coastal  management  strategies. 

Quality  of  life 

The  improvements  to  coastal  wave  prediction  models  as  developed  in  this  project,  will  improve 
modeling  capability  of  coastal  circulation,  morphological  development,  surge  prediction  and  transport 
processes,  which  will  benefit  coastal  recreation  (more  reliable  knowledge  of  wave  heights,  rip  currents 
etc),  coastal  management,  and  help  mitigate  pollution  hazards  for  humans  (recreation)  and  coastal 
ecosystems. 

TRANSITIONS 

Economic  Development 

The  developments  in  this  project  will  be  made  available  as  open  source  software  and  as  modules  to 
widely  used  operational  wave  models.  These  models  are  used  by  NOAA/NCEP  and  other  agencies 
involved  in  coastal  development  and  management,  and  by  many  coastal  and  ocean  engineering 
companies. 

Quality  of  life 

The  software  developed  within  this  project  will  be  disseminated  in  open  source  models  used  by  local 
and  federal  agencies  and  companies  involved  in  coastal  recreation  (surf  prediction,  rip  currents, 
pollution,  surge  prediction),  coastal  management,  and  mitigation  of  coastal  hazards. 

RELATED  PROJECTS 

Coastal  Wave  Observations  at  FRF,  Kitty  Hawk,  NC,  USA.  This  project  is  carried  out  by  Jeff  Hanson, 
Kent,  Hathaway  and  Harry  Friebel.  It  is  strongly  related  to  our  project  for  exchange  of  field  data. 

Modeling  Wind  Wave  Evolution  from  Deep  to  Shallow  Water;  Nonlineary  and  Dissipation.  Grant 
N0014-10-1-0453.  Pi’s:  Tim  Janssen  (San  Erancisco  State  University),  Tom  Herbers  (Naval 
Postgraduate  School)  and  Gerbrant  van  Vledder  (Delft  University  of  Technology). 

SWAN  and  SWASH  development  teams.  Delft  University  of  Technology,  Delft,  The  Netherlands. 
http://www.swan.tudelft.nl  &  http://swash.sourceforge.net. 
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